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The spectrum of meson and diquark excitations in cold color-superconducting (2SC) quark matter 
is investigated under local color and electric neutrality constraints with /3-equilibrium. A 2-flavored 
Nambu-Jona-Lasinio type model including a baryon /j,b, color [is, and electric /j,q chemical poten- 
tials is used. The contribution from free electrons to the free energy is added to take into account the 
/3-equilibrium. The sensitivity of the model to the tuning of the interaction constants in the diquark 
(H) and quark-antiquark (G) channels is examined for two different parameterization schemes by 
choosing the ratio H/G to be 3/4 and 1, respectively. At H — 3G/4 the gapless- and at H = G the 
gapped neutral color superconductivity is realized. It is shown that color and electrical neutrality 
together with /3-equilibrium lead to a strong mass splitting within the pion isotriplet in the 2SC 
phase (both gapped and gapless) , in contrast with non-neutral matter. The n- and a- meson masses 
are evaluated to be ~ 300 MeV. It is also shown that the properties of the physical S'[/(2) c -singlet 
diquark excitation in the 2SC ground state varies for different parameterization schemes. Thus, for 
' H — 3G/4 one finds a heavy resonance with mass ~ 1100 MeV in the non-neutral (gapped) case, 

whereas, if neutrality is imposed, a stable diquark with mass ~ \^q\ ~ 200 MeV appears in the 
gapless 2SC environment. For a stronger attraction in the diquark channel (H — G), there is again 
a resonance (with the mass ~ 300 MeV) in the neutral gapped 2SC phase. Hence, the existence 
of the stable massive SU(2) c -singlet diquark excitation is a new peculiarity of the gapless 2SC. In 
addition, the behaviour of the diquark mass in vacuum, i. e. , at fiB — 0, as a function of H has 
been investigated. 
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O ■ I. INTRODUCTION 

According to modern theoretical observations made in the framework of perturbative QCD, at asymptotically 
high baryonic densities and low temperatures the strongly interacting quark matter is expected to undergo a phase 
transition to the color superconducting state [j], Q ■ Unfortunately, a perturbative QCD analysis is not applicable at 
moderate baryon densities (which might exist inside compact stars or in heavy ion collision experiments) and a study 
is usually done with the help of effective theories, such as the Nambu-Jona-Lasinio (N JL) model [1, 0, Q . 

In spite of the lack of quark confinement in NJL, it successfully describes low-energy pseudoscalar and vector mesons 
^ . in the hadronic phase (see, e. g., @,0|)- This success is provided by the fact that many of light meson properties, e. g., 
masses, are driven by chiral symmetry, rather than by confinement. Moreover, the chiral phase transition at high 
temperatures and/or density expected in QCD, is naturally described by the NJL model 0, H|. The NJL model is 
also well-suited for the consideration of a hot and/or dense medium under the influence of external conditions [9. Hof 
as well as for the investigation of different physical processes in it [ll|, OjJ . 

In the earlier studies of color superconductivity [HI, [3, [H| for the case of two-flavor quark matter (u and d quarks), 
only the influence of baryonic density was taken into account. From these investigations, it became evident that the 
two- flavor color superconducting phase (2SC) might be yet present at rather small values of \xb ~ 1 GeV, i. e., at 
baryon densities only several times larger than the density of ordinary nuclear matter (see reviews 

0QEH2I)- This 

is just the density of compact star cores. 

The quark matter inside compact stars is considered to be electrically and color neutral in a bulk. Moreover, 
there must be an equilibrium between the gain and loss in the /3-decay d — ► u + e + v e (here e is the electron 
and v e is the electron neutrino), the so-called /3-equilibrium. All these physical constraints must be taken into 
account when studying the equation of states for a compact star. To do this in the NJL model, additional chemical 
potentials related to the electric charge density of quarks and electrons as well as to color charges must be introduced. 
Recently, an intensive theoretical study of neutral color superconducting quark matter has been given (see, e. g., 



Kecently, an intensive tneoretical study 

[H [H El [H, HJ HI HI HI IH [H HI), which revealed a new possible ground state of the 2SC phase, where some 
additional number of quasiparticles with a gapless dispersion law is appeared [l!| (it is the so-called gapless color 
superconductivity g2SC, an antipode to the usual, gapped 2SC). 

There is a great interest in the study of different excitations of the 2SC phase, as the application of its results 
in some related investigations, e. g. in astrophysics and heavy-ion collision experiments it may help to reveal some 
observable effects evidencing the formation of a quark-gluon plasma. In particular, the bosonic excitations of the 2SC 
phase ground state, such as tt- and er-mesons as well as diquarks, are expected to be copiously produced in dense 
medium with rather strong correlations between quarks and antiquarks (just this thing is realized at moderate baryonic 
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densities) and affect some scattering and decay processes to a visible effect. Diquarks on their own are also of great 
interest in hadron physics because of their importance in determining baryon properties (27j . As to the compact stars, 
one may find an influence of these particles on the equation of state and on the cooling process. Moreover, in dense 
matter, the deconfinement phase transition might be accompanied by the appearance of Bose-Einstein condensed 
diquark matter [28[ etc. 

In our recent papers [IS OS], we have studied the masses of mesons and diquarks that are formed in cold (T = 0) 
and dense quark matter in the framework of a 2-flavored NJL model with baryon chemical potential fis- In particular, 
it was shown that in the 2SC phase the meson masses lie in the interval 330-^500 MeV, depending on the values of 
fiB G (1050, 1200) MeV. Since the original SU(3) C color symmetry of the model is spontaneously broken down in 
this phase to SU(2) C , one may expect the appearance of five Nambu-Goldstone bosons. However, we have proved 
that the abnormal number of three, instead of five, massless bosons is allowed for the diquark sector of the model. 
In addition, there are two light diquarks as well as a heavy diquark resonance that is an SU(2) c -singlet with the 
mass ~ 1100 MeV. Qualitatively, the local color neutrality constraint does not affect the masses of mesons and the 
SU(2) c -singlet diquark. However, in this case only one Nambu-Goldstone boson and four light diquarks are present 
in the 2SC phase in the NJL model [3l[ (see also the discussion at the end of the paper). At nonzero temperature 
some of the properties of mesons and diquaks in a strongly interacting quark matter were discussed in [32l |33| . 

In the present paper, we continue our investigation of mass spectra for mesons and diquarks, imposing the local 
electrical neutrality and /3-equilibrium to the cold 2SC medium, in addition to the color neutrality. As in our previous 
papers, we use a 2-flavored NJL, where additional color (/i 8 ) and electric ((J-q) chemical potentials are introduced. It 
will be shown that in the 2SC phase, the color neutrality constraint supplemented by the electrical neutrality and 
/3-equilibrium drastically changes the mass spectrum of the 7T-, a- and diquark, in comparison with non-neutral quark 
matter. 

The paper is organized as follows. In Section II, the thermodynamic potential as well as the effective action of the 
NJL model, extended with baryon (/is), color (us), and electric (/xq) chemical potentials, are obtained in the one-loop 
approximation in /3-equilibrium. Further, in Section III, the gap equations and the phase diagram of quark matter 
are investigated under the local color and electrical neutrality constraints. Here, the behaviour of fi$ and /iq vs \xb 
are obtained for neutral quark matter with 2SC type color superconductivity within two different parameterization 
schemes: H — 3G/4 and H = G. In the first case, H — 3G/4, a gapless 2SC phase revealed itself, whereas for H = G 
a gapped phase is preferred. In Sections IV and V, some peculiarities of the mass spectra of the it, <j mesons and 
scalar diquarks are investigated both in the gapless and gapped neutral 2SC phases. (In addition, the influence of the 
diquark channel coupling constant on the diquark mass in the vacuum, i. e. at [Lb = 0, is also considered.) Finally, 
in the Appendix, the expression for the quark propagator in the Nambu-Gorkov representation is obtained. 



II. THE MODEL AND THE EFFECTIVE ACTION 



Our investigation is based on the NJL type model with two quark flavors. Its Lagrangian describes the interaction 
in the quark-antiquark as well as scalar diquark channels: 



L„ 



'fidv — m 



G 



{qqf + {qil 5 fq) 2 \+H ^ [fif 5 T 2 X A q] [qi~f 5 T 2 \ A q C } , 
,4=2,5,7 



(1) 



where the quark field q = qi a is a flavor doublet (i = 1,2 or i = u,d) and color triplet (a = 1,2,3 or a = r,g,b) 
as well as a four-component Dirac spinor; q c = Of and cf = q l C are charge-conjugated spinors, and C = ij 2 j° 
is the charge conjugation matrix (the symbol t denotes the transposition operation). It is supposed that up and 
down quarks have an equal current (bare) mass m. Furthermore, r a stands for Pauli matrices, and for Gell- 
Mann matrices in flavor and color space, respectively. Clearly, the Lagrangian L q is invariant under transformations 
from color SU(3) C as well as baryon U(l)s groups. In addition, at m = this Lagrangian is invariant under the 
chiral SU(2)i,xSU(2) j r group. At m ^ the chiral symmetry is broken to the diagonal isospin subgroup SU(2)/ 
with the generators Ik = T&/2 (fc = 1,2,3). Moreover, in our system the electric charge is conserved, too, since 
Q = I3 + B/2, where I3 is the third generator of the isospin group SU(2)j, Q is the electric charge generator, and B 
is the baryon charge generator (evidently, these quantities are unit matrices in color space, but in flavor space they 
are Q — diag(2/3, — 1/3), -Z3 = diag(l/2,—l/2) and B — dia<?(l/3, 1/3)). If the Lagrangian dXJ) is obtained from the 
QCD one-gluon exchange approximation, then H = 3G/4. In addition, we find it interesting to deal with another 
relation between coupling constants, H — G, which results in qualitatively different model properties (see below). 

In order to take into account /3-equilibrium, we include electrons into our consideration, extending the Lagrangian 
as follows 



L 



qe 



(2) 



Here e is the electron spinor field (for simplicity, electrons are taken to be massless). Clearly, the Lagrangian ([2]) is 
well-suited for the description of different processes in the vacuum, i. e. in the empty space. Since the prime object of 
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the present paper is the consideration of dense medium properties, we extend the Lagrangian ((2|) by including terms 
with charge densities and chemical potentials as it is usually done in statistical physics 

L = L qe + n B N B + n Q N Q +fi 8 N 8 . (3) 

In (|3|), Nb, Nq, Ns are baryon, electric and 8i/i-color charge density expressions, respectively; /is, fJ-Q, fJ-8 are the 
corresponding chemical potentials. Recall that 

N B = qB~f°q, Nq = q&fq - e 7 °e, N s = qT 3l °q, (4) 
where T 8 = \/3A 8 — diag(l, 1, —2) is a matrix in the color space. From equations ([¥]), we have 

VbN b + HqNq + fj, 8 N 8 = H e ej°e + ^ M M <fe7 C Vfa = M e e7°e + ?M7% (5) 

where fi e is the electron number chemical potential, and fj,i a is the chemical potential for the number of quarks with 
color a and flavor i. Obviously, one has 

_ Ms 2^i Q A*b MQ i 

M ur — M«3 — — H g r , Mdr — Mdg = "~g" 7^ + M8 j 

Mb 2^q Ms MQ 
Mub = — + -g 2^i 8 , /ijj = — — -2fig, fj, e = -fi Q , (6) 

A = Mb/3 + ^qQ + ^ 8 T 8 = M + <fytr 3 + /i 8 T 8 , (7) 

where the last equality in |[7J) is obtained due to the above mentioned relation Q = I^+B/2; moreover M = 

5fx = /xq/2. It follows from (J6j) that [i^a — M«a + Me for each color a. The matrix M is diagonal in the six-dimensional 

(color) x (flavor) space, and its matrix elements are just the quantities Hi a from (|6j). 

If all chemical potentials in ([3]) are nonzero and independent quantities, then SU(3) C and SU(2)j are not the 
symmetry groups of this Lagrangian. Instead, due to the /i 8 - and ^(Q-terms, it is symmetric under the reduced 
color SU(2) c xU(1)a s and flavor U(l)j 3 groups. With the local neutrality imposed, the chemical potentials in Q 
are, however, no more independent quantities of the model. Equating further (Nq) and (Ng) to zero, the chemical 
potentials are subjected to two constraints, thereby fixing two of them. As a result, fj,Q and /x 8 become dependent on 
fiB- It turns out (see below) that there exists a critical value of the baryon chemical potential [i B in the locally neutral 
matter which separates two phases: if /j,b < Ms the normal quark matter phase with /xq = and fig — is formed, 
and the Lagrangian ([3]) is an SU(3) c xSU(2)/ invariant one; if [Ib > Hb (2SC phase), the fiQ and ^ 8 are not already 
equal to zero and have a nontrivial /is-dependence. It means that at \ib — (J-b the SU(3) c xSU(2)j symmetry of the 
Lagrangian (J3]) is explicitly (not spontaneously) broken by the chemical potential terms (i. e. no Nambu-Goldstone 
bosons must appear) to the color SU(2) c xU(1)a 8 and flavor U(l)j 3 groups. (Note that the color U(1)a s group is 
broken spontaneously in the 2SC phase, see below.) Due to the above-mentioned flavor symmetry transformation in 
the critical point n B , one could expect that all pions would have equal masses for /i_b < /i^, whereas at larger \xb 
(fiB > (J-b) the pion mass splitting should occur in neutral matter. 

To study the phase diagram of the system and the mass spectra of meson and diquark excitations, we need to get 
the thermodynamic potential as well as an effective action up to second order for the bosonic degrees of freedom. Since 
electrons and quarks are not mixing, the total thermodynamic potential f2 of the system is the sum of its electronic tt e 
and quark fl q parts: fl = fl q + fl e . It is well-known that fl e — —[i\jY2'K 1 . To obtain f2 g , we start from the Lagrangian 
describing the quark contribution only (see ([2])), 

C = L q + qfrfq, (8) 

where fi is the quark number chemical potential matrix, defined in ([7]) . The linearized version of Lagrangian (jHJ) that 
contains auxiliary bosonic fields has the following form 



C = q 7 I/ i<9 !/ + M7° - cr-m - ij 5 Tr a T a 



a a + 7r a 7T„ 



- ±A* A A A -^[q C il 5 T 2 \ A q}-^[q il 5 T 2 \ Aq C ], (9) 

where here and later, a summation over repeated indices a — 1,2,3 and A, A' — 2,5,7 is implied. Clearly, the 
Lagrangians (|SJ) and ([9]) are equivalent, as can be seen by using the equations of motion for bosonic fields, which take 
the form 

a(x) = -2G(qq), n a (x) = -2G(^7 5 t q <z), A a (x) = -2H(q c rf^X^), A* A {x) = -2H{ qil 5 T 2 X Aq c ). (10) 

One can easily see from ()10|) that the mesonic fields <j(x), ir a (x) are real quantities, i. e. (a(x)y = cr(x), (7r a (x)y — 
ir a (x) (the superscript symbol f denotes the hermitian conjugation), but all diquark fields A^(x) are complex scalars, 
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so (Aa(x))^ = A* A (x). Clearly, the real cr(x) and ir a (x) fields are color singlets, whereas scalar diquarks Aa(x) form a 
color antitriplet 3 C of the SU(3) C group. If some of the scalar diquark fields have a nonzero ground state expectation 
value, i. e. (A^a;)) ^ 0, the color symmetry of the model (|8|) is spontaneously broken down. 

It is more convenient to perform our investigations in terms of the semi-bosonized Lagrangian @ , since in this case 
we have a common footing for obtaining both the thermodynamic potential and the effective action of the model. 
Indeed, in the one fermion-loop approximation, the effective action S c s(a, 7r a , Aa, A a , ) of the model (J9j) is expressed 
by means of the path integral over quark fields: 



exp(iS cS (a, n a , A a, A* a ,)) = N' J [dq] [dq] exp (i J td 4 



where 



S c b(u, 7T a , A^, A^,) 



d 4 x 



4G 



A A 

AH 



S, 



off ■ 



(11) 



and N' is a normalization constant. The quark contribution to the effective action, i. e. the term S c s in (JTTJ) , is given 
by: 



exp(iS cff ) = N' J [dq][dq] exp(^- 
In (jT2J) we have used the following notations 



qD + q + q^D-q u - qKq c - q° K* q 



d 4 x 



(12) 



D+ = vfd v - m + #7° - S, D~ =i-i v d v -m-lii ) -Y}, £ = a(x) + i~/ 5 w a (x)T a , 
= a (x) + i 1 5 i: a {x)T t a , K* = iA^(x) 7 5 T 2 A^, K = ^A J 4(x)7 5 T 2 A J 4, 



(13) 



where are nontrivial operators in coordinate, spinor, color and flavor spaces. 1 In the following, it is very 
convenient to use the Nambu-Gorkov formalism, in which quarks are composed into a bispinor W such that 



* = 



, & = y = (q,q C ) = (q,q t C) = y t 



C , 1 



Furthermore, by introducing the matrix-valued operator 

Z = 



D+, -K 
-K*, D~ 



(14) 



(15) 



one can rewrite the functional gaussian integral in (|12p in terms of \t and Z and then evaluate it as follows (clearly, 
in this case [dg][dg] = [dq c ][dq] = [d*]): 



exp(i<S off ) = / [d*] exp 



*Z^d 4 ; 



= / [d*] exp \ - I *'(FZ)$d 4 .x \ = det 1/2 (rZ) = dct 1/2 (Z), 



where the last equality is valid due to the evident relation detY" = 1. Then, using the general formula detO = 
expTrlnO, one obtains the expression for the effective action: 



<Seff(c, "a) A A, A* A ,) 



d 4 X 



&aA* a 



4G 



AH 



-Tr 



sfcxNG 



In Z. 



(16) 



Besides of an evident trace over the two-dimensional Nambu - Gorkov (NG) matrix, the trace in (JTHJ) is calculated in 
spinor (s), flavor (/), color (c) and four-dimensional coordinate (x) spaces, respectively. 

Starting from (|16[) . one can define the quark contribution Q q (cr, 7r a , A^, A^,) to the thermodynamic potential (TDP) 
of the model (|8|). In the mean-field approximation one has 



cr,7r a , ,A* , —const 



-n q (a, n a ,A A ,A%) J d 4 



(17) 



1 In order to bring the quark sector of the Lagrangian (0 to the expression, given in the square brackets of 1121 . we use the following 
well-known relations: d* = -d v , C^C- 1 = -(7")*, C^C~ Y = (7 s ) 4 = f 5 , t 2 tt 2 = -(f)*, r 2 = ^ °. ' ~M. 
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The ground state expectation values (mean values) of the fields: (a a (x)) = a°, (ir a (x)) = tt°, (Aa(x)) = 
A A , (A^, (x)) = A A °, are solutions of the gap equations for the TDP fl q (in our approach all ground state ex- 
pectation values do not depend on coordinates x): 



dn 1 = Q 

dlta 



da 



o. 



dn q 

MI 



0, 



dn q 



0. 



(18) 



Next, let us perform the following shift of bosonic fields in ([TBI : <j{x) — > a(x) + er°, Tr a (x) — » ir a (x) + tt°, A a (x) — > 
A^(x) + A^°, A A (x) — > A j4 (a;) + AJ4. (Obviously, the new shifted bosonic fields a(x),ir a (x), A A (x), A^(x) now 
denote the small quantum fluctuations around the mean values a°, 7r°, A° A , A* A ° of mesons and diquarks rather than 
the original fields (fTU)) ). In this case 



Z = 



-K„ 



-K* D~ 



E , K 
K* , E* 



E , K 
K* , E* 



(19) 



where So is the quark propagator matrix in the Nambu - Gorkov representation (its matrix elements Sij are given in 
the Appendix IBj) . and 

(K , K*,Df, E , E*) = (tf, £>± E, E ( ) 



(7 — CT . 7T a — 7T„ 



Then, expanding the obtained expression into a Taylor-series up to second order of small bosonic fluctuations, we 
have 



S eB (a,7r a ,A A ,A A ,) = 5$ +S^(a,7r a ,A A ,A A ,) + ---, 
where (due to the gap equations, the linear term in meson and diquark fields is absent in l|20p) 



: (2) 



(20) 



S, 



(o) 
cir 



(fx 



A A A* A ° 



S^(a,7r a ,A Al A* A ,) 



d 4 x 



<7 2 + 7T, 



ag ah 

= -n q (<T°,ir° a ,A A ,A* A °,) J fx 
2 A A A^ 



^ scfxNG ln(S x ) 



AG 



:Tr 



AH 



scfxNG 



S 



E 
K* 



K 

E* 



So 



E 
K* 



K 

E' 



(21) 



(22) 



In the following we will study the spectrum of meson/diquark excitations in the color superconducting phase of the 

(2) _____ 

N JL model on the basis of the effective action S\J . The effective action (|2"2"|) can be presented in the explicit form: 



where 





c(2) _ 
^eff - 


= 

^mesons 




c(2) 
diquarks 




c(2) 
mixed 


-^ v scfx {5'llE5i2 



5(2) +s m +S {2) 

-'mesons 1 ^diquarks 1 mixed' 



-Tt scfx {S 12 K*S 12 K* + 2S 11 KS 22 K* + S 21 KS 21 K} . 



(23) 

(24) 
(25) 
(26) 



and Sij are the matrix elements of the quark propagator matrix So defined in (|19p (see also Appendix |BJ). Moreover, 
some necessary explanations concerning the trace-operation over coordinate space in the expressions (|24[) - (|26 () are 
given in Appendix [A] (see (|A4jl ). It follows from these formulae that the effective action (f24| is a functional of the 
meson fields <?(x), iti{x) only, the effective action (|25f is composed from diquark fields only, and the mixing between 
mesons and diquarks might occur because of the effective action (f2"6"| . 



III. GAP EQUATIONS AND NEUTRALITY CONDITIONS 

Earlier (see, e.g., the papers (l9l. l20l. I2H] ) . it was shown that for electrically and color neutral cold dense matter, 
described in the framework of the model ([3]), only two phases are allowed to exist. In the first one, the symmetric 
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phase that is usually called the normal quark matter phase, only the mean value of the c-field, (c(x)) = a°, is 
nonzero. In the second one which is just the 2SC phase of dense matter the mean value of the diquark field A 2 (x) 
is nonzero in addition, i. e. (A2(x)) = A 2 7^ 0. Hence, without loss of generality, it is convenient to deal with TDP 
S}q(a°, 7r°, A° A: A*fi) (note, A, A' — 2,5,7), in which all arguments, except a° = M — m (m is a bare quark mass, 
whereas the parameter M is usually called constituent or dynamical quark mass) and A 2 = A, are identically equal to 
zero. In this case the calculation of the total TDP £1 of the system ([3]) (recall, Q, = Vt e + O g ) is significantly simplified 
fl9L [20II , and we have 



, (M-mf 



12tt 2 AG AH 

- 2 E / 4^{\ E i + 5 ^ + \ E i-^\}-T. J ^{\E±^u b \ + \E±f, db \}, (27) 



where E± = ^/{E± f + |A| 2 , £± = E ± /2, E = V^ + M 2 , /I = ( Mur + /i dg )/2 =(/x„ s + /i dr )/2 =/x + /x 8 =/x B /3 + 
Hq/6 + /i 8 , /i„b = + 2^q/3 - 2^i 8 , ^ df) = /is/3 - ^q/3 - 2^ 8 , <5^ = /iq/2 (see also the notations, used in (|7|)). 

Note that, apart from the order parameters M and A, the TDP f2 really depends on the chemical potentials, which 
is indicated in |27|) in an explicit form. Since the integrals in the right hand side of (|27|) are ultraviolet divergent, 
we regularize them as well as the other three-dimensional divergent integrals below by implementing a cutoff in the 
integration regions, |q| < A. Starting from (j27|) . one can find the gap equations dil/dA* — and dfl/dAI = 0, which 
supply us with the values of M, A in the ground state of the system: 



OA* 
dM 



AH 
M — m 
2G 



2A 



d 3 q 
(2tt) 3 



6(E+ - |M) , 0(E A - \6fM\) 



E 



A 



2M 
AM 



d 3 q 
(2ir) 3 E 

d 3 p 
(2irfE 



0{E- \ f i ub \) + 9(E- \fx db \) 
8(E+-\5n\)E+ , 9(E- £ 



El 



+ 



0. 



\6n\)E- 



E\ 



= 0. 



(28) 



(29) 



Next, let us impose the local color- as well as electric charge neutrality requirements on the ground state of the model 
(|3|). It means that the quantities fj,s and [iq take such values that the densities of the 8-th color charge N$ and 
electric charge Nq are equal to zero in the ground state for arbitrary fixed values of other model parameters, i. e. 
(N$) = —dfl/dfis = 0, (Nq) = —dfl/dfiQ = 0. These neutrality constraints look like: 



d 3 c 



e(E+ 



M)E+ 9(E A - M)E- 



-4 



(Nq) 



3tt 5 



(2tt) 3 
d 3 p 
(2tt) 3 

(N s 



El 



E\ 



+ 



sign(n ub )0(\[i ub \ - E) + sign(n db )9(\fj, db \ - E) 
d 3 q 



= 0, 



6 



sign(fi ub )6(\fi ub \ - E) 



sign(6n)9(\5fi\ - El) + sign(5fi)8(\5fi\ - E A ) = 0, 



(30) 



(31) 



In the following we suppose, for simplicity, that A is a real quantity. Of course, the solution of the common system of 
gap equations (|2"51) - (|2"9")l and neutrality relations (f5U|) - (f5T)) is possible only numerically. In all numerical calculations 
of the present paper we use the following parameter set 



G = 5.86 GeV~ 



A = 618 MeV, 



= 5.67 MeV 



(32) 



that leads in the framework of the NJL model to the well-known vacuum phenomenological values of the pion weak- 
decay constant = 92.4 MeV, pion mass M w = 140 MeV, and chiral quark condensate (qq) — —(245 MeV) 3 . 
Moreover, we use two different values for the coupling constant H in the diquark channel, H = 3G/4 and H = G, for 
which two qualitatively different 2SC phases are realized in the model (see below) . The numerical analysis shows that 
for the parameter set (|32[) and both relations H = 3G/A and H = G the system of equations (|2"8"|) -([3"T j) has only two 
solutions. As was already discussed after ([7]), the first one (with M ^ 0, A = 0, = and = 0) corresponds to 
the SU(3) c xSU(2)/ symmetry of the model (normal phase), whereas the second one (with M / 0, A / 0, /i8 / and 
fiQ 0) corresponds to the 2SC phase. As usual, solutions of these equations give local extrema of the thermodynamic 
potential fi(/iB) A*q, A*s; M, A) (|27p. Clearly, one should also check in which of them the TDP takes the least value 
(for each fixed value of /ib)- Just this solution of the system of equations (f2"B")) - (|3Tj) corresponds to the genuine neutral 
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FIG. 1: The behaviour of M vs /i_g in the neutral matter. 



FIG. 2: The behaviour of A vs \ib in the neutral matter. 



ground state of the model. In particular, for each fixed lib it supplies us with mean values (gaps) M, A (see Figs 1,2) 
as well as with values /Zg and liq (see Figs 3,4), at which the ground state has zero charges. 2 

It is clear from Figs 1,2 that there is a critical value of the baryon chemical potential \i c B (for the case H = 3G/4 
one has ll c b — 1.08 GeV, whereas for the case H — G one gets \i c B = 1.04 GeV 3 ) such that at /ib < ^°b the normal 
SU(3) c xSU(2)/-symmetric phase of the model occurs. However, at hb > Ms the 2SC phase is realized. Since in this 
case liq 7^ (see Fig. 4), the isospin SU(2)/-symmetry of the normal phase is broken in the critical point rather 
by hand than dynamically, down to the U(l)/ 3 group. Furthermore, in the 2SC phase we have (A 2 (a;)) = A ^ 0, 
(Ac,j(x)) = 0. So in the critical point ll c b the color SU(3) c -symmetry is broken down to the SU(2) C subgroup. Since 
at [Lb > /J-b the /zs-term of the Lagrangian Q is nonzero (see Fig. 3), the initial color symmetry of the model is also 
broken by hand down to the SU(2) c xU(1)a s subgroup, which is further broken dynamically (spontaneously) down to 
SU(2) C . Hence, one may expect the appearance of only one Nambu - Goldstone boson in the mass spectrum of the 
2SC phase (see the discussion below). Finally, remark that in the critical point fi B the transition between these two 
phases is of the first-order which is characterized by a discontinuity in the behavior of M and A vs. lib (see Figs 1,2). 

Up to now we have discussed the properties of the 2SC phase in neutral and /3-equilibrated matter which are common 
for the two particular cases H = 3G/4 and H = G. However, there exist, at least two, qualitative distinctions between 
these neutral 2SC phases. (Both of them are based on the fact that in the neutral 2SC matter the relation A < \Sfi\ is 
valid for the case of H = 3G/4, whereas at H = G the opposite one A > |<5^| is true.) The first one lies in the diquark 
mass spectrum and will be discussed below. Just now we would like to present the second one which is provided 
by the quasiparticle dispersion relations, i. e. the momentum dependence of energy. In condensed matter physics 
quasiparticles are simply the one-fermion excitations of the ground state. In our case the quasiparticle spectrum of 
the 2SC matter is defined by singularities of the quark propagator Sq (see Appendix [Bj) . Clearly, there are twelve 
(6 quark and 6 antiquark) quasiparticles in the 2SC matter. Four of them, blue quasiparticles, have the energies 
E ± Li u b and E ± /idb- The energy spectrum of the other eight, red and green quasiparticles, consists of four values 
E^ ± \Sfx\, each is doubly degenerate. Evidently, for both relations between coupling constants H and G there are 
momentum values at which the energies of two blue quasiparticles must turn into zero, i. e. there are no energy costs 
to create these quasiparticles (their energies are E — Li u b and E — /i c ib)- Due to this reason, these excitations are called 
gapless ones. Similarly, at H — 3G/4 there are in addition also two gapless red and green fermionic excitations (with 
energies E^ — |<5/i|) of neutral 2SC matter, since in this case the relation A < \S/i\ is true. However, in the neutral 
2SC matter of the case H = G there are no additional gapless red and green quasiparticles. Hence, at H = 3G/4 the 
quasiparticle spectrum of the neutral 2SC matter consists of four gapless excitations (it is a so-called gapless color 
superconductivity). In contrast, at H = G there are only two gapless quasiparticles in the neutral 2SC matter, and 
color superconductivity is called as gapped one. Thus, our results corroborate the conclusion, made, e.g., in [19l.l23l]. 
that the gapless 2SC may exist for a rather narrow interval of the coupling constant H values. 

In the next sections we will calculate the inverse two-point (unnormalized) correlators of meson and diquark fluctu- 



2 In the literature, the zero value of the current quark mass m is of frequent use in the 2SC investigation. In this case the constituent 
quark mass M is identically equal to zero in the 2SC phase of neutral matter (see, e.g., in [Tsl |22| V 

3 The tendency, the gre ater H the smaller fig , is indeed supported by earlier investigations of the 2SC phenomenon in the framework 
of NJL models |34l, |35|I . In particular, it was shown that at sufficiently large values of the coupling constant H the 2SC phase may be 
realized in the model even at zero baryon chemical potential iig |35| . 
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FIG. 3: The behaviour of us vs [ib in the neutral matter. 



FIG. 4: The behaviour of \iq vs \xb in the neutral matter. 



ations over the ground state of neutral 2SC matter in the one-loop (mean-field) approximation and find their masses 
in the two cases, H = 3G/4 and H = G. 



IV. MESON MASSES 



After a more detailed study of the S^ xcd -part (|26|) of the effective action, it turns out that it is composed from 
cr(x), A 2 (a;) and A^cc) fields only, i. e. the 7r-mesons are not mixed with diquarks (this property is justified by the 
parity conservation both in the normal- and 2SC phases of our model) . So to find the 7r-meson masses it is enough to 

deal with the effective action <Smesons (l24l) . which is a generating functional of the one-particle irreducible (IPI) Green 
functions of a- and 7r-meson fields. In this case, instead of 7r i (cc)-fields ([TU]) . we will use the new fields ir°(x) = tts(x), 
7r ± (a;) = (7Ti(a;) ± 7T2(x))/v2, so that the quantities E, E* from (I13|) look like: 



<j{x) + ij 5 {TT°(x)T3 + 7T (x)t + + TT + (x)T-\ 



a(x) + Z7 5 [7r° (x)t^ + ir (x)r_ + 7r + (x)r+] 



(33) 



where t± = (ri ± r 2 )/y/2. Then, the IPI Green functions of a- and n^, 7r°-meson fields can be generated through the 
relation: 



XY 



(x-y) 



w Ernes' 



SY(y)5X(x) ' 



(34) 



where A, F = cr, , 7r° (to take the variational derivatives in ([34]) , it is very instructive to refer to the relations 
(IA4[) - (IA5|l ). In momentum space the zeros of the Fourier transformations of these functions are connected with meson 
masses. 



A. The 7r ± -meson masses 

Using the relation ({M)) . it is possible to define the IPI Green functions of the 7r ± -fields. In particular, 

IW(a; -V) = ^ + 2 * Tr .sc [2A 11 (z) 7 5 S 11 (- 2 ) 7 5 + 2B 22 (z) 7 5 A 22 (-z) 7 5 
+A 12 {z) 1 b A 21 {~z) 1 5 + B 21 (z) 1 5 B 12 {-z) 1 5 ] 



(35) 



(in this expression all traces over flavor indices are calculated, and z — x — y), where operators Aij (z) and Bki(z) are 
defined in Appendix IB1 Moreover, T 7T - 7T + [x — y) = T^+^- (y — x) and T 7T + 7T + (x — y) = r^-^- (x — y) = 0. Due to the 

traces, containing an odd number of 7 , all the mixed Green functions of the form r cr7r o,± (x — y) — — STri>.±(y)8™(x) are 
zero. Moreover, since the traces containing an odd number of r±-matrices equal zero, all the IPI Green functions of 
the form T^o^i (x — y) are also zero. Hence, in the framework of our model there is no any mixing between Tr^-fields 
on the one hand, and ct, 7r°-fields on the other. After taking the trace over color indices, the Fourier transform of (|35p 
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has the following form: 



1 



r\r+,r-(p) = +4iTr s 



2bg 2) (P + 9)7 5 «^ J (?)7 5 + (P + 9)7 5 «r 2 J (9)7 5 + 3S(p + g) 7 W(g)7 5 + 621 (p + <zh 5 M<z)7 5 



,,5.(12), 



(2tt)* 
>,(3) ( 



2a£ a) (p + «)7 5 ^ 1 i 2) (9)7 5 + ag>(p + ?)'M! ) (?)7 B + 



„5„(3), 



(36) 



where the quantities a l:L (q), b^ (q) etc are the corresponding Fourier transformations of a\i (z), (z) etc, 
presented in Appendix [Bl Clearly, T w - w +(p) =T w + w -(—p). The zeros of these functions determine the 7r ± -meson 
dispersion laws, i. e. the relations between their energy and three-momenta. In the present paper, we are mainly 
interested in the investigation of the modification of meson and diquark masses in dense and cold color- and electrically 
neutral matter. Since in this case a particle mass is defined as the value of its energy in the rest frame, p = (see, 
e.g., (33, HH, we put p = (poj 0,0,0) in the following. As a result, the calculation of 1PI Green functions is 

significantly simplified. Indeed, in the rest frame one can easily perform all the trace calculations over spinor indices 
in (IBU)) (see the auxiliary relations IjBlOp ) and gets 



1 



IVtt- (Po) = ^ - 16« 



A (q -S}i)(pa+q +5n)-E + E- - |A| 



-4i 



d A q 



(2tt) 4 {( qo - 5fif 

1 



(£ A ) 2 ][(p + q + 5fl)* 



(El) 2 } 



1 



(Po + qo + M«6 + E)(q + fi d b - E) (po + qo + M«6 - E)(q + fi db + E) 



(37) 



where we have used the same notations as in ([27]) . Note also that in ([37|) qo is a shorthand notation for qo + ie -sign (go) 
and (po + go) is a shorthand notation for (po + go) + • sign(po + go) , where e — > 0+ (see also the remark after (|B9[) and 
[13] )• The go-integration in (f37|) is performed along the real axis in the complex go-plane. We will close this contour 
by an infinite arc in the upper half of the complex go-plane. Inside the obtained closed contour the integrand of the 
first integral in (|37p has four simple poles which are located in the following points: 



(g )i = 5fi - E A + is ■ 0(E A - 5fi), 
(g ) 3 = E+ - 5fi - po + is ■ 0(5/1 - £+), 



(90)2 = 5fi + E A +ie- 0(-E A - 5fi), 
(qo)i = -E^ -5(i-p + ie- 9(5(1 + El), 

whereas the integrand in the second line of ([37]) has the following four poles in the upper half of the go-plane: 

(g )i = E - Hdb + i£ ■ 6((i db - E), (g ) 2 
(90)3 = ~E - fi u b - Po + is ■ 0(E + fiub), (qo)i 

Summing the residues of the integrand function in these poles, we can perform the g -integration in (|37[) and obtain: 



-fidb - E + ie ■ 6(fMb + E), 
E - fiub -po+ie ■ 6 (fi ub - E). 



rV +7r -(po) = 
6(E A -5fi) (po-E A 



d 3 



q 



(2tt) 3 
25fi)E A + 



(-5fi -E A ) (p 



E 



A 



2Sfi)E l 



E+E- 



IAI 



E 



A 



(po + E A + 25fi) 2 - (E+) 



+ \2 



E+E- 



IAI 



6(5(i 



El) 



(Po 



El 



25(i)El 



E+E- 



E\ 



(po -E A + 25(iY - (Elf 



El 



6(5(i + El) ( Po + E+ + 25fi)E+- E+E- 



El 



(po 



-El 
6(E 



25 fi) 2 
fidb) 



(e- a ? 

6(fl u b 



(p Q -El + 25fif-(E A Y 
6(fi db - E) 



d 3 g 

(27T) 



Po + 25fi + 2E 



E) 



6(fi ub + E) 

pa + 25(i + 2E ^ p + 25fi-2E p n + 25fi-2E 



(38) 



Clearly, the matrix element T„+ n - (po) depends effectively on the variable z = (po + 25 fi). It follows from our 
numerical analysis that for both values of the coupling constant H (H — G and H — 3G/4) the expression (j3"5|) 
has only two zeros, zi(fis) and z 2 (M-b)- Hence, for each fixed value of \ib > 11%, i- e. in the 2SC phase, we have 
r7r+7r-(Po) ~ (po + 25fi — zi)(po + 25fi — Z2). Since (po) = r w + T -(— po), the determinant of the inverse propagator 

matrix of the 7r ± -mesons has the following form: 



r w + w -(po) ■ r 7r - 7r +(p ) ~ (p - (25(i - zi) )(p - (25(1 - z 2 ) ). 



(39) 



Evidently, in the pp-plane it turns into zero in two points, which are the mass squared of the 7r ± -mesons. Hence, in 
the 2SC neutral dense matter the 7r ± -mesons have different masses: 



Ml + = (25(i - Zl f, Ml. = (25fi - z 2 f 



(40) 



The behaviour of M w + and M n - vs fiB in the color superconducting and neutral matter is depicted in Figs 5,6 for 
the cases H = 3G/4 and H = G, respectively. 



10 




FIG. 5: The behaviour of meson masses vs \ib in the 
gapless (H = 3G/4) color superconducting neutral matter 
(Ha > Mb = 1-08 GeV). 



FIG. 6: The behaviour of meson masses vs fis in the 
gapped (H = G) color superconducting neutral matter 
Ob > Mb = 1-04 GeV). 



B. The n , (7-meson masses 



As it was already discussed after ([55)1 . 7r -mesons are not mixed with other fields in the framework of the NJL 
model ©, i- e. the 1PI Green functions of the form T^xix — y), where X{x) = a(x),7T ± (x), Aa(x), A\,(x), are 
equal to zero. In this case the 1PI Green function T 7T a 7: o(x — y) is just the inverse 7r°-propagator, which can be found 
from (|3"4")l . To get the 7r°-mass we need the expression r 7r o 7r o(p), i. e. the Fourier transform of r„.o 7r o(a; — y), in the 
rest frame. Using the technique presented in the previous section, one can obtain after tedious but straightforward 
calculations at p = (po,0, 0, 0): 



1 



I\r°7r°(Po) - 2G 

d 3 q fe(5n-E A 



- 8 



d 3 q E+E A + E+E- + A 2 E+ + E A 



(2tt)< 



(2tt) 3 

-5fj,-E A ] 



E A E A 



E 



A 



pI-(e+ + e a y 

(po + E A )E A -E+E- -|A| 2 (£ A -p )£ A -£ + £- 



(Pc 



^2 



6(5fJ, - E£) + 6(-5ii - E} 



-8 



A 

E 



(E+ - PQ )E+-E+E-- 
(E + A -p f-{E A f 



IA 



12 f 17+ I ..\Pt 7^+ TTn- I A 121 



2 - (^) 2 



(E++ PQ )E+-E+E- -|A| 2 



rf 3 g 

(2tt) 3 p 2 - 4S 2 



(i?++Po) 2 -(^ A ) 2 
-/i db - E)} . 



(41) 



The function (|4T|) is an even one, i. e. r„.o„.o(po) — ^tt°tt° ( — Po)- Hence, it effectively depends on the variable p\. 
Numerical investigations, performed both at H = 3G/4 and H = G, show that for each fixed value of /is > Ms the 
1PI Green function (|4"Tj) has a single zero on the positive p 2 ,- semi- axis, which is just the mass squared of the 7r°-meson. 
Its mass in the 2SC phase of the dense and neutral matter is presented graphically in Figs 5,6. 

The situation with the c-meson mass is much more involved. Indeed, as it follows from the previous section, in 
the model under consideration the a- meson is mixed with A2(x), Aj (a;)-diquarks (see also [2{| [3(1 HH). So to get 
the particle masses in this case, one should find the zeros of the 3x3-matrix determinant, whose matrix elements are 
nothing but the 1PI Green functions Txy(po) (we use the rest frame in the momentum space representation), where 
X(x),Y(x) = <t(x), A2(x), A|(:e). In general, it is a rather hard task, which however can be significantly simplified 
due to some reasons. It turns out that 1PI Green functions of the form T a x (X(x) = A2(x), A^x)) are proportional 
to the gap A (see Fig. 2) and constituent quark mass M (see Fig.l) as well. Hence, in the normal color symmetric 
phase (where A = 0) there is no mixing between the a-meson and diquarks Aa^), A^x) at all. The similar is true for 
the 2SC phase, if the current quark mass m is zero (in this case, as was pointed out in the footnote^ the parameter M 
is also equal to zero in the 2SC phase) . In our consideration the parameter M is a rather small quantity in comparison 
to the gap A in the 2SC phase (see Figs 1,2). So to have a grasp of the order of magnitude of the er-meson mass, we 
ignore for simplicity the a — A2 mixing effect in this case, too. As a result, we see that both for H — 3G/4 and H = G 
the cr-meson mass is defined by the 1PI Green function r CTCr (a; — y) which is approximately the 1PI Green function 
r^o^o [x — y). Hence, for both above mentioned values of the coupling constant H the cr-meson mass is approximately 
equal to the 7r°-meson mass in the 2SC neutral matter (see Figs 5,6). 
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Now some comments, concerning the behaviour of the 7r-meson masses in the neutral dense matter, are in order. 
First of all, note that in the normal phase (where fig < ^ C B ) the /^Q-term of the Lagrangian ([3]) is zero. As a result, 
the ground state of this phase is an SU(2)/-invariant one. Due to this symmetry, in the normal phase all 7r-mesons 
have a common mass that is approximately 140 MeV for all hb < Hb both for H — 3G/4 and H = G (the a-meson 
mass is approximately equal to 700 MeV in the normal phase). Above the critical point, i. e. at [Ib > Mbj the isotopic 
SU(2)j symmetry of the system is broken due to the appearance in (|3|) of a nonzero /ig-term. So in the 2SC phase of 
electrically neutral matter the meson masses are allowed to have different values. Just this conclusion was supported 
by our numerical investigations (see Figs. 5,6), where a rather strong splitting of 7r-meson masses is observed. In 
contrast, if the electric charge neutrality is not imposed on the NJL-system, then all mesons have a common mass in 
the 2SC phase [H[3l|. 

Finally, note that in [38[ a rather strong splitting of 7r-meson masses was shown to exist in electrically neutral and 
noncolor-superconducting dense quark matter both with or without pion condensation phenomenon. 



V. DIQUARK MASSES 



As in the previous section, we will ignore, for simplicity, the mixture between diquarks A2(x), Ajj^x) and er-meson 
(this is justified by moderately small values of the parameter M in the 2SC phase). In this case, in order to obtain 

( 2^ 

the masses of diquarks, we need to analyze the 1PI Green functions, generated by the effective action Sd iC j Uarks (|25|) : 

r "<- 9) =-Iot- (42) 

where X(x),Y(x) = A A (x) , A* A , (x) . 



A. Diquark masses in the 2SC phase 

Due to the structure of (|25p . the diquarks, as such, are not mixed to one another in the framework of our model. 
So it is reasonable to study step-by-step the diquark excitations of the 2SC neutral matter in the A2(x), Aj-^x)-, 
A 5 (x), Ag(x)-, and finally in the Ay(x), A^x) sectors of the model. 

The investigations of the 1PI Green functions (|42|) in the A$(x) and At(x) sectors of the model supply us with four 
excitations of the 2SC phase ground state. All of them have the common mass 3|/xs| both for H = 3G/4 and H = G, 
when color- and electric charge neutrality constraints are imposed [22|, |31| . Evidently, these excitations form two real 
(or one complex) doublets of the SU(2) C ground state symmetry group of the 2SC phase. 

To study the masses of the SU(2) c -singlet diquark excitations, we should consider the 1PI Green functions in the 
A2(x), A^(x) sector of the model. It can be shown in the usual way that in the 2SC phase these quantities take the 
following form (here again the rest frame, i. e. p = 0, in the momentum space representation is used): 



r A2 A 2 (po) = r A .A 2 -(po) = 4A 2 / (Po)> r A .A 2 (po) = r A2 A*(-po) = (4A 2 - 2^)I (^) + Ap h{pi), (43) 
where 



Wo) 



,,_ / d 3 q 9(E+ - |^|) f d 3 q e(E^-\Sfi\) 
(27r)3£+ [4(25+)2-p§] J (2tt)3 £-[4(^)2 -pj 

, , d 3 q 9{E+ - \5y\)E+ f d 3 q 0(E A - \S^l\)E- 



d 3 g 9(E+ - \Sfi\)E+ f d 3 g 9(E A - |fy|)£T 
(2ir) 3 E+[4(E+f-p>] ' ( 27 3 E~ A [A{E A Y - pi] 



From the expressions (|43[) it is possible to compose the inverse propagator matrix Q 1 (po) for the diquarks 
A2(x),A|(x) moving in neutral 2SC matter: 



0-1 (p„) = - 1 JWpo). r^Cpo) j (45) 



r A * A2 (po), r A | A; (p ) 

Then, the mass spectrum in the A2 (x) , A2 (x) sector of the model is defined by the equation 

deta- 1 ^) = *pl{(pl - 4 A 2 )/o(Po) - ^ 2 M)} = ^PlF(pl) = 0. (46) 



In the pp-plane this equation has an evident zero, corresponding to a Nambu-Goldstone boson (NG), p\ = 0. (Since 
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FIG. 7: The sketch of the SU(2) c -singlet diquark res- 
onance in the gapped (H = G) 2SC neutral matter 
(/is > A*s — 1-04 GeV). The solid line is for its mass Mo 
vs hb, whereas the width of the strip between dashed 
lines Md ± T/2 is its width r vs /j,b- 
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FIG. 8: The behaviour of the stable SU(2) c -singlet di- 
quark mass Md vs /j,b in the gap less (H — 3G/4) color 
superconducting neutral matter (fiB > A*s = 1-08 GeV). 



F(0) ^ 0, it does not have another zeros of the form p 2 , = 0.) 4 Since in the 2SC phase the chemical potential //g is 
nonzero, it is clear that in this phase the initial color SU(2) c xU(1)a 8 symmetry of the Lagrangian ([3]) is spontaneously 
broken down to SU(2) C group. So, there is only one broken symmetry generator, corresponding to the above mentioned 
NG-boson solution of the equation (|4"5)) . Hence, it is possible to assert that in the 2SC phase of color- and electrically 
neutral matter, described by the NJL model ([3]), there appears a normal number of NG-bosons, corresponding to 
a spontaneous breaking of the color symmetry. In contrast, if neutrality requirements are not imposed and only the 
baryon chemical potential is taken into account, then in the NJL model there is an abnormal number of NG-bosons 
in the mass spectrum of the 2SC phase 29]. (Note, an abnormal number of NG-bosons is not a quite unexpected 
phenomenon. It is inherent to a variety of quantum models with broken Lorentz symmetry, which is provided by 
chemical potentials [28l 1391].) 

For the further investigation of the equation (|4"rj)l some additional information about the functions Io,Ii (|44j) . and 
consequently about the function F in (|46[) . is required. It turns out that these functions are analytical in the complex 
Pg-plane, except for the cut n < p\ along the real axis (the first Riemann sheet). It is easy to verify that je = 4A 2 
in the case of the gapped 2SC phase, i. e. at H = G, where A > \6fi\. However, x = 4|<5p| 2 = fig > 4A 2 in the 
case of the gapless 2SC phase, i. e. at H = 3G/4, where A < \8p\. The quantity F(pq) is a complex- valued function 
in the whole first Riemann sheet, except for points on the real axis which do not belong to a cut, where F(j)q) is a 
real-valued function. 

Notice that in the case of the gapless 2SC phase, the cut of the pg-plane originates to the right of the point 
4A 2 , whereas in the case of the gapped 2SC phase it just starts in the point 4A 2 . This circumstance is of decisive 
importance for the appearance of the diquark mass difference in the two types of color superconductivity. Indeed, at 
H = G we did not manage to find any solution of the equation F(pq) = in the first Riemann sheet of the variable 
Pq. (In particular, it is evident from ([46]) that in the real points such that pg < 4A 2 the function F(pf ) ) is an exactly 
negative quantity.) So using the procedure, presented in the appendix of our previous paper [3(|, we continue the 
function F(pq) to the second Riemann sheet, where it takes a zero value in some complex point. It means that an 
SU(2) c -singlet diquark resonance appeared in the mass spectrum of the neutral gapped 2SC phase. Its mass Mo and 
width r are presented in Fig. 7 as functions of /is. For the gapless 2SC phase (at H = 3G/4) the situation is quite 
different. In this case the first Riemann sheet of the function F{pf ) ) contains in addition the set of real points pg such 
that 4A 2 < Pg < Pq. Just among these points the zero M|> of F(pq) is located. It means that the existence of an 
SU(2) c -singlet stable diquark excitation with mass Md such that 2A < Mo < ImqI ( see Fig- 8) is typical for the 
neutral gapless 2SC phase (at H = 3G/4), in contrast to the neutral gapped one (at H — G), where it is a resonance 
(see Fig. 7). In addition, it is easily concluded that for the same relation H = 3G/4 the SU(2) c -singlet diquark is a 
heavy resonance in the 2SC (gapped) phase without neutrality requirement (30l . l3ll | , whereas it is a stable particle in 



4 It might seem that the appearance of the NG-boson in the mass spectrum of the model is strongly connected with the disregarding of 
the mixing between cr-meson and A2(x), AS (a;) diquarks. However, as it was shown in [30ll . this NG-boson is a native property of the 
model, since it exists in the mass spectrum even in the case, when the mixing is taken into account. 
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the gapless neutral 2SC phase, if the neutrality constraints are fulfilled. 5 

Notice that in the present paper actually the dispersion relations of mesons and diquarks, i. e. the momentum 
dependence of their e nerg ies, are investigated for the particular value of the three-momentum, p = 0. Contrary, in 
the recent papers [24], |25|] the diquark dispersion relations were studied at p ^ 0, and a conclusion about the Higgs 
instability of the NJL gapless 2SC phase was made. In this case the g2SC phase is unstable vs fluctuations of the 
Nambu - Goldstone diquark fields and, as a result some of the diquarks acquire a non-physical negative velocity 
squared (24|. (In addition, the g'2SC phase of a gauged NJL model suffers also from a chromomagnetic instability 
[261].) In particular, it was shown that a sufficient condition for the Higgs instability in the gapped 2SC phase is the 
relation v2|^mI > A. Since in the case G = H this relation is not fulfilled, one may conclude that in the gapped 
neutral 2SC phase considered above the Higgs instability is absent. 

Finally, it is necessary to note that in the physical SU(2) c -singlet diquark channel, we have found a stable massive 
mode in the case of the g2SC phase (H = 3G/4), whereas in 24] it was claimed that a gapless tachyon emerges in 
this channel. Actually, there is no contradiction between our result and that paper. Indeed, in our case the diquark 
mass Mjj obeys the relation 2 A < Mb = po- In contrast, in [24| those solutions for the diquark dispersion relation 
were investigated that are constrained by IpoMpI ^ A. It means that in the gapless color superconductor there 
are two branches for the dispersion relation of the physical SU(2) c -singlet diquark mode at small values of \p\. The 
first is the massive heavy excitation with 2 A < po (see Fig. 8 of the present paper), the other one corresponds to a 
gapless tachyon, located in a quite different kinematic region with |po| A. In addition, it was shown in [2~ij that 
the singularity po/\p\ appears in the diquark two-point 1PI Green function in the gapless 2SC phase. In this case, 
owing to this singularity, it is quite possible to get at \p\ — > a result (gapless tachyon) that does not coincide with 
the diquark mass, directly calculated at p = 0. 



B. Diquarks in the normal phase (A = 0, [i$ = 0, /iq = 0) 

At A = the three complex diquark fields Aa(x) (A = 2, 5, 7) are not mixed with other fields in the second order 
effective action (|23l) of the model. Since in addition the quantities fig, /iq are equal to zero in the normal phase, there 
is an SU(3) C symmetry of the ground state of this phase. So, in order to study the diquark masses at fis < Ms; ft 
is enough to consider, e.g., the A2-diquark sector only. In this phase the inverse propagator matrix G~ 1 (j>o) for the 
diquarks A2(x), A^(x) looks, in contrast to the inverse propagator (|43|) of these fields in the 2SC phase, more simpler 
[3(j,[3l| (again, we use the rest frame, p = 0, in the momentum space representation): 



- 1 (po) = -( 



o , r A2 A*(po) 



rA*A 2 (po), o 



(47) 



where r A2 A*(po) = r A*A 2 (-po), 



= W- 16 J Wf 4Ei-(p E + 2» B /^ S IE - *«' (48) 

and e = (po + 2/xb/3) 2 . Since the determinant of the inverse propagator matrix (|47|) takes the form 

dct^Opo) = T^~ 2 (p )T^ 2 (p ) = T^~ 2 (po)TZ~Z~ 2 {-p Q ), (49) 

it is clear that the mass spectrum is defined by the zeros of the 1PI Green function (148|) . Note, the function <E>(e) is 
analytical in the whole complex e-plane, except for the cut AM 2 < e along the real axis. (In general, this function is 
defined on a complex Riemann surface which is to be described by several sheets. The integral representation for $(e), 
given in (j48| . defines its values on the first sheet only. To find a value of $(e) on the rest of the Riemann surface, a 
special procedure of analytical continuation is needed (see, e.g., in 30]).) It turns out that for the model parameter 
set (I3"2l and a wide set of the coupling constant H values (see below), the equation 

*(e) - ^ (50) 



5 Let us quote also another less rigorous, but more physical argument in favour of the stability of the SU(2) c -singlet diquark excitation 
in the g2SC phase (H = 3G/4). Namely, the decay of the above mentioned diquark excitation with mass Mp into a pair of u- and d 
quark quasiparticles is forbidden due to the energy disbalance. Indeed, the total energy E(p) of this quark pair in their center-of-mass 
system takes the form E(p) = E u (p) + Ed(p), where p = \p\ > and E u (p) = E^ + \Sfi\, E^{p) = |E7 — \Sfi\\. One can easily show that 
E{p) > ImqI- Since the mass Mo of the SU(2) c -singlet diquark excitation in the g2SC phase is smaller than |jUq|, we conclude that the 
decay of this excitation into a pair of u- and d quasiparticles is forbidden. 
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FIG. 9: The behavior of the stable diquark mass M£> vs 77 = H/G in the vacuum (at [ib = 0). Here 77* = H* /G ~ 0.698, 
77** = H** /G ~ 1.525, M ~ 0.350 GeV. At 77 < rf there are no stable diquarks in the normal phase, and these particles are 
resonances. At 77 > 77** an SU(3) C symmetric ground state of the normal phase (including the vacuum) is unstable in favor of 
color superconducting phases. 



has on the first Riemann sheet of the variable e a single root eo on the real axis such that < eo < AM 2 . In this case 
the physical meaning of eo is that eo = (Af£,) 2 , where is the mass of a stable diquark at /is = 0. Having a root 
eo, one can find two zeros of the 1PI Green function Ta*a 2 (po) as well as four zeros of the determinant (|4"5|) . They 
provide us with the following mass squared in the AJj, A2 sector of the model: 

(Ma) 2 = {M° D - 2 MB /3) 2 , (M A .) 2 = {M° D + 2,i B /3) 2 . (51) 

In particular, if H = 3G/4, then M° D w 1.988M, if H = G, then M° D « 1.746M (here M is the constituent quark 
mass, or gap, in the normal phase (see Fig. 1), i. e. M « 0.350 GeV). We relate Ma in (f5"Tj) to the mass of the diquark 
with the baryon number B = 2/3 and Ma* to the mass of the antidiquark with B = —2/3. The difference between 
diquark and antidiquark masses in (|51[) is explained by the absence of a charge conjugation symmetry in the presence 
of a chemical potential fis- 

Finally, due to the underlying color SU(3) C symmetry, the previous statement is valid also for Ag, A5 and Ay, A7. 
As a result, we have a color antitriplet of diquarks with the mass Ma (l5Tj) as well as a color triplet of antidiquarks 
with the mass A/a* in the mass spectrum of the normal phase, i. e. at [Ib < 

It is clear from the equation ([50]) that its own solution eo lies inside the interval < eo < 4M 2 only if H * < H < H** , 
where H* and H** are defined by 

1 

H* = 



4$(4Af 2 ) 4 [AV M 2 + A 2 + M 2 ln((A + VM 2 + A 2 ) /M)] ' 

4$(0) 4 [AVM 2 + A 2 - M 2 ln((A + VAf 2 + A 2 )/Af)] 2(M - m) 1 ' 

(here the values of A, m and G are presented in ([32]) . whereas M w 0.350 GeV is the dynamical quark mass in the 
normal phase). In this case, as was noted above, eo defines the stable diquark mass M£, in the vacuum, i. e. at [ib = 0, 
through the relation eo = (Af£>) 2 . The values of M£, are depicted in Fig. 9 as a function of the parameter 77 = H/G. 
(Note, 77* = H*/G, 1]** = H**/G =1.5M/(M - m).) For a rather weak interaction in the diquark channel (H < H* 
or 77 < 77*), eo runs onto the second Riemann sheet, and unstable diquark modes (resonances) appear. Unlike this, a 
sufficiently strong interaction in the diquark channel (H > H** or 77 > 77**) pushes eo towards the negative semi-axis, 
i. e. (Af£,) 2 < in this case. The latter indicates a tachyon singularity in the diquark propagator, evidencing that the 
SU(3) c -color symmetric ground state is not stable. Indeed, at a very large H, as it has been shown in [35| . the color 
symmetry is spontaneously broken even at a vanishing chemical potential. The fact that at 77 — > 77" the diquark mass 
Mjj tends to zero may be considered as a precursor of the spontaneous breaking of the SU(3) C symmetry, taking part 
a±H = H**. 



VI. SUMMARY AND DISCUSSION 



The present paper is the last in the series of our papers [29J, |30j, 13_1| , devoted to the investigation of the ground state 
bosonic excitations (mesons and diquarks) of color superconducting quark matter. The novel features of our present 
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consideration, performed in the framework of the two-flavored NJL model, are the local electric charge neutrality 
constraint as well as the /^-equilibrium that is due to taking into account electrons. As a result, depending on the 
diquark channel coupling constant H, the gapless (H = 3G/4) or gapped (H = G) 2SC may exist in the system 
based on the Lagrangian ([3]). Since in ([3]) a new term with the electric charge chemical potential /iq appeared, we 
have an explicit breaking of the flavor SU(2) symmetry in both 2SC phases. So, in contrast to the non-neutral case 
[lij [111 [HI] , 7r-mesons acquire a rather strong mass splitting in both gapless- and gapped 2SC phases (see Figs. 5,6). 

The diquark sector of the model consists of six modes. Since in the color and electrically neutral 2SC phase the 
chemical potential ji% is not equal to zero (see Fig. 3), the Lagrangian §5§ is invariant under the color SU(2) c xU(1)a 8 
group. It turns out that in the color superconducting phase this symmetry is spontaneously broken down to SU(2) C , 
so, in accordance with general theorems, one of the diquark modes takes zero mass, i. e. it is the Nambu - Goldstone 
boson that is an SU(2) c -singlet. Besides, in this phase there are four very light diquark excitations that are composed 
into two real SU(2) c -doublets with a common mass proportional to fig. In our opinion, when the NJL model is gauged, 
the above mentioned diquark modes should be absent in the model. Instead, due to the Anderson - Higgs mechanism, 
five massive gluons will appear. So we believe that five diquark excitations with zero or very small mass, observed in 
the 2SC phase of the nongauged model ([3]), are not physically interesting objects and might be ignored (see also the 
discussion in [32J). 

The remaining one, an SU(2) c -singlet diquark excitation of the 2SC ground state, has different properties under 
different external conditions. Indeed, if the electric charge neutrality requirement is not imposed, then, as it was 
shown in [30 . [31| at H — 3G/4, this diquark mode is a heavy resonance with mass approximately equal to 1100 
MeV. However, if the neutrality and /3-equilibrium conditions are imposed, then, on the one hand, at H = 3G/4 the 
gapless 2SC phase is realized, in which the above mentioned diquark mode is already a stable particle, whose mass is 
evaluated around 200 MeV (see Fig. 8). Note that its mass is approximately equal to the value of |/xq| presented by 
Fig. 4. On the other hand, at H — G we have the neutral gapped 2SC phase, in which the diquark is still a resonance 
but with a much more smaller mass (see Fig. 7). 

It was shown, e.g., in [40j that at sufficiently high baryon densities, comparable with densities inside neutron stars, 
the normal quark matter equation of state is significantly influenced by scalar mesons (er-meson etc). Moreover, these 
mesons might change significantly the mass and radius of the neutron star as well as the role of hyperon degrees 
of freedom in dense matter. Since the masses of diquarks in the electrically neutral color superconducting quark 
matter are of the same order as those of mesons (compare Figs 5,6 and Figs 7,8), one might expect similar effects in 
the measurable parameters of neutron stars (e.g. their masses, radii etc), when diquarks will be taken into account 
properly in neutron star physics (see also [28jj). So, the influence of rather light diquark excitations might in principle 
be checked by astrophysical observations. 

In fact, in the present paper the dispersion relations for mesons and diquarks following from the NJL model © are 
investigated in the rest frame, i.e. at \p\ = 0. In the recent paper [24| (see also (25|) the diquark dispersion relations 
at small nonzero \p\ were studied in the same model. It was found there that in the gapless 2SC phase (H — 3G/4) 
Nambu - Goldstone diquark modes have a non-physical negative velocity squared, i.e. the model is unstable against 
fluctuations of these fields. (In contrast, as it follows from our discussion at the end of section IV A[ in the gapped 
2SC phase (H = G) this kind of instability is absent.) Moreover, it was proved in [24[ that at small |p| ^ in the 
SU(2) c -singlet diquark channel of the g2SC phase there is a gapless tachyon with |po| <C A. Our result, i. e. the 
existence of a stable diquark excitation with mass Mo ~ 200 MeV (see Fig. 8) in the same channel, does not conflict 
with this, since po = Mo belongs to a quite another kinematic region, where po = Mo > 2A. In spite of the fact that 
in our calculations \p\ — 0, we believe that at \p\ ^ there is also a branch of the SU(2) c -singlet diquark dispersion 
relation, which corresponds to a stable excitation of the g2SC phase with mass Mo (see Fig. 8). 

Moreover, in the present paper the behaviour of the diquark mass at vanishing /j,b as a function of the coupling 
constant H is also obtained (see Fig. 9). 
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APPENDIX A: SOME FORMULAE 

Appendix A contains some useful formulae employed in the text, 
i) Determinant: 



d.'i ( £ ' J ) = det l~ CB + CAC^D] = det[DA - DBD~ l C\. (Al) 



ii) Inverse matrix: 



A, B\ 1 _ / C~ X DL , -N \ _ f L ,_ -A^BN 
C , D I ~ [ -L , B~ l AN I ~ I -D~ X CL , N 



(A2) 



where 

i = [y4C _1 i) - B] -1 , N =[DB- 1 A-C}- 1 , Z = [A- BD^C]' 1 , AT = [D - CA^B]' 1 . (A3) 

iii) Variational derivatives: Let A, B are some operators in the coordinate space with matrix elements 
A(x,y) =A(x — y) and B(x,y) = B{x — y), respectively. Moreover, let <j(x) and 4>{x) are some fields. Then, 



Tr{AaB(j)} = J dxdydzduA(x, z)cr(z)8(z — y)B(y,u)(j)(u)8(u — x) = j dxdyA(x,y)cr(y)B(y,x)<f>(x). (A4) 
It follows from fA4)) that 

6 2 Tr{AaB(j)} 



8<j{y)5cj){x) 



A(x,y)B(y,x) = A(x - y)B(y - x). (A5) 



APPENDIX B: QUARK PROPAGATOR IN THE NAMBU-GORKOV REPRESENTATION 

In the Nambu - Gorkov representation the inverse quark propa gat or matrix Sq 1 is defined in (|19[) . Using the 
relation (|A2|) as well as the energy projection operator technique of [34[, one can obtain the following expressions for 
the matrix elements Sij(x — y) of the quark propagator So(z) (here z = x — y): 

Sn(z) = An(z)r + r_ + Bh(z)t_t + ; Si 2 (z) = A 12 (z)t + + B 12 (z)t_; 

S 22 {z) = A 22 {z)t + t_ + B 22 {z)t_t +] S 21 (z) = A 21 (z)r + + B 21 {z)t_, (Bl) 

where their explicit structure in the flavor space is presented with the help of t± = (t\ ± t 2 )/^/2 matrices (recall, t% 
are the Pauli matrices) and 

All(Z) ^i(2^ e 1 (g + W -(^) 2 7 A+ + (^W 7 / " 

2 7 (2tt) 4 \g +£ + M«& qo + Vub-E) 3 11 U 12 11 W 3 ' V ; 

+ I / /^e-^ ( 7 ° A+ + 7 ° A - 1 P 3 ^ ee b[\ 2 \z)pg + b[l\z)Pt\ (B3) 



i 4 ^ 

2 7 (2tt)4 c \ (g,, _ 5 M )2 _ {Ea) 2 ' J1 + ^ (eo _ 5/i) 2 _ (jB + )2 

1 f dS _ lQZ f 7 °A+ , 7 °A- 



2 7 (2tt) 4 



{ 7 \ = + 7 ^ J ^ - og») WPfiJ + «g> WP^, (B4) 
L QO — Hub + E go - Hub — & J 



5 22 (z) = I [j£^ e -*>{ (m + W-E- - (q + Sh) + E+ . (c) 

^ 22(Z) 2 7 (Jhr)*' 1 (g + ff M )a - (i^ 7 A+ + (g + Stf ~ (EX)^ 1 1 



2 
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AM = /If L + + WTW^ } s a '* ,Aj ' (B6) 



</2 J OY I (40 - Jfi) 2 - (E^f (,„ - 4>) 2 - (Elf 

*■« " ^/(&-'" L + ,ffW + " < B8 > 

, , , A*A 2 f d 4 q inr f 7 5 A+ 7 5 A_ ] . . 

A 21 (z) = =i / <^ /.„ + - ... + — — - \ = a 21 (z)X 2 . B9 

y/2 J (2tt) 4 I (q - fy) 2 - {E+f (qo — Sfi) 2 — (E A ) 2 J 

In the above formulae P[f =diag(l, 1, 0), P% =diag(0, 0, 1) are the projectors on the red-green and blue subspaces of 
the color space, correspondingly; A2 is the Gell-Mann matrix; A± = 7 S£ij£) ) are projectors on the solutions of 

the Dirac equation with positive/negative energy. The other notations appearing in (|B2j) - (|B9|) are identical to those 
of (|27p . Note that in (|B2)) - (|B9)) and similar integrals, containing an integration over the energy variable, the symbol 
go is a shorthand notation for qo + is ■ sign(go)) where e — ► 0+. This prescription correctly implements the roles of 
fiB, Ms and fJ-Q as chemical potentials and preserves the causality of the theory (see, e.g. |37|). Introducing the new 

projectors A± = i(l ± Xil&Ml^ it is very convenient to use in trace calculations the following relations 

7 5 A±7 5 = A±, 7°A±7° = A T , A 2 ± = A±, A±A T = 0, 

TrA ± -2, Tr(A ± A ± ) = ^, Tr(A±A T ) = —j^-- (BIO) 
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